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Abstract: 



^ \ We show that near a point where the equihbrium density of eigenvalues of a matrix 

model behaves like y ~ x^^'', the correlation functions of a random matrix, are, to 
leading order in the appropriate scaling, given by determinants of the universal 
^ . {p, g)-minimal models kernels. Those (p, q) kernels are written in terms of functions 

^ ■ solutions of a linear equation of order q, with polynomial coefficients of degree < p. 

^ ■ For example, near a regular edge y ~ x^^"^, the (1, 2) kernel is the Airy kernel. Those 

^ I kernels are associated to the (p, q) minimal model, i.e. the [p, q) reduction of the KP 

O \ hierarchy solution of the string equation. Here we consider only the 1-matrix model, 

O ■ ior which q = 2. 

> 

X 

Contents 

1 Introduction 

1.0.1 Example: (1,2) law: Airy kernel and Tracy- Widom law 3 

1.0.2 Higher (p, 2) laws 14 

1.1 Universality of eigenvalues statistics point of view . . . 

1.1.1 Double scaling limit 

1.1.2 Specific heat 

1.2 Formal matrix models and combinatorics point of view 



^ E-mail: michel.bergere@cea.fr 
^ E-mail: bertrand.eynard@cea.fr 



1 



1.2.1 Spectral invariants and spectral curve 

1.2.2 Limits of large maps 

1.2.3 Singularities of spectral invariants 

2 Minimal models 

2.1 String equation 

2.2 Tau function 

2.3 Lax pair 

2.4 Lax equation 

2.5 The differential system 

2.6 Correlators 

2.7 Example: Airy kernel 

2.8 Classical limit 

2.9 Topological expansion 

2.10 BKW expansion 

2.11 Symplectic invariants 

2.12 Double scaling limit and (p, 2) kernel 

2.13 Parametrics of orthogonal polynomials and Baker- Akhiezer functions 

3 Kontsevich's integral 

4 Derivatives 

4.0.1 Example 

4.0.2 Example W^^'^ 

4.0.3 Other examples 

4.0.4 Homogeneity relation 

5 The (2m, 1) minimal model 

5.1 The (2m, 1) minimal model 

6 Conclusion 



1 Introduction 

In this article, we shall consider "scaling limits" of matrix integrals. 

We shall show, under certain assumptions, that scaling limits of matrix integrals are 
governed by some well known integrable systems. The fact that double scaling limits of 
matrix models are minimal models (p, q) of conformal field theories [26], has been well 
known in the physics literature for a long time (see [27l [20t |4T] for review, and among 
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others see [5ll IMl SHI SIl [53 [52 [IS [23 |13 [23l |12] ) here we merely summarize 
some results scattered in the physics literature, we present the main features of those 
universal limit laws, and provide a mathematical proof. 

The idea of the proof works backwards: we show that {p, q) minimal models deter- 
minantal correlation functions satisfy the same recursion as the scaling limits of matrix 
models. 

We shall consider only the 1-matrix model, whose corresponding limit integrable 
systems are the (p, 2) minimal models, reductions of KdV, and we hope to later gener- 
alize those results to multi-matrix models and general (p, q) limits, as claimed in many 
physics works [271li8l[5T] . 



The main result, theorem 12.31 is that limit correlation functions are given by deter- 
minantal formulae of the (p, q) kernel. 

1.0.1 Example: (1,2) law: Airy kernel and Tracy- Widom law 

The equilibrium density of eigenvalues of a x random hermitian matrix, generically 
behaves near the edge of the distribution, like: 

p{x) ~ x^. (1-1) 

It is well known that, after rescaling x by A^^''^, the ra— points correlation functions in 
the vicinity of the edge, are given by determinants of the Airy kernel which appears in 
Tracy- Widom law [66] of extreme eigenvalues statistics: 



p„(A-2/3xi,...,A-2/3a;„) ~ A^det(i^Airy(x„x,)) (l + 0(A-^/3))^ (1_2) 

N—*oo 



and the Airy kernel is the ChristofFel-Darboux kernel of the Airy function: 

AAiry(Xi,X2) = . (1-3) 

X\ — X2 

Notice that the Airy function satisfies a 2nd order ODE, whose coefficients are poly- 
nomials of degree 1: 

Ai"{x)=xAi{x), (1-4) 
which can also be written as a 2 x 2 differential system: 

d_ ( Ai{x)\^ [0 1\ f At{x)\ 

dx [At'ix) J [x Oj [At'ix) J ■ ^''^^ 
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Figure 1: Generically, the equilibrium density of eigenvalues behaves like x^/^ near 
endpoints. It may happen, after fine tuning the parameters, that it behaves like x^^'^ 
{q = 2 for 1-matrix model). The eigenvalue statistics at a scale x = x N'^/^'p^'^^ is 
governed by the universal (p, q) kernel law. The (p, q) kernel is an integrable kernel, 
associated to the (p, q) reduction of the KP hierarchy. For a regular endpoint (p, q) = 
(1, 2) it is the Airy kernel. 



1.0.2 Higher (p, 2) laws 

More generally consider a p/q singularity of the equilibrium density of eigenvalues: 

p{x) ~ xP/". (1-6) 

We shall consider only q = 2 and p = 2m + 1 in this article, but we recall that the 
physics literature claims that general y ~ x^^*^ case can be treated the same way, and 
should correspond to (p, q) minimal models. 

We shall see, after rescaling x by A^'J/lp+g)^ ^j-^^t the correlation functions in the 
vicinity of the edge, are given by determinants of the (p, q) kernel which appears in 
(p, q) minimal models of conformal field theory 



Pn{N-l^.xi,...,N-l^^Xn) ~ iV"^det(i^(p,,)(a:„x,)) (l + 0(iV^i/(P+«))), (1-7) 



and the (p, 2) kernel is the Christoffel-Darboux kernel of the (p, 2) Baker- Akhiezer 
function: 

'llj{Xi)ij{x2) -1p{Xi)ij{x2) ,^ 



K(p,2){Xl,X2) 



Xi - X2 
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where the [ip, tp) functions satisfy a 2nd order ODE: 

where T>^^p2){x) is a 2 x 2 matrix with polynomial coefficients, such that the degree of 
det V is of degree at most p. 

Moreover this differential system, is associated to the Lax matrix of the {p, 2) reduc- 
tion of the integrable KdV hierarchy, which means that the coefficients of the matrix 
are themselves solution of some non-linear integrable differential equations. 

The coefficients of T'(p_g)(a;) are differential polynomials of a function u{t) which 
satisfies (for p = 2m + 1, g = 2), the m + 1^^ Gelfand-Dikii non linear equation [27], of 
the form: 

The time t measures the distance to the critical point. 

For example for pure gravity {p,q) = (3,2), V(^3 2){x) is a 2 x 2 matrix, with poly- 
nomial coefficients such that the degree of det P(3,2)(a^) is 3: 

^(3'^)(^'') = ((x + 2n)(2x-2«)+n ^"^'t) ^^"^^^ 
and where u{t) is solution of the Painleve I equation (1*^* Gelfand-Dikii equation): 

3u^it)-^ = t. (1-12) 

All this has been stated for a long time in the physics literature, and we shall just 
present it concisely and prove it. 

1.1 Universality of eigenvalues statistics point of view 

In this subsection, we summarize some well known facts about random matrices [571 
EOl E31 El ESI EZl [131 [m , and we fix the notations. 

Consider a probability law of the form of the ioint law of eigenvalues of a random 
hemitia„typ<imatr,x: 

dA<(Ai,...,A„) = i J](A, -A,)= Y[e-T'-'i>:> d\u (ngj 

i<j i 



■^Hermitian matrices correspond to real eigenvalues and positive measure d/i, but it is customary 
to generalize random matrices to normal matrices having their eigenvalues on some contours in the 
complex plane, and the measure dfi can be complex. The loop equations are the same for all those 
models, they are independent of the integration contour, and thus, they can all be treated in the same 
framework. 
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where Z is the partition function 

Z = ^ dfi= jdM e-^*=-^W. 



:i-i4) 



Here s is a parameter, often called the temperature, we shall be interested in the large 
N limit, and possibly a limit s — > Sc, where Z has a singularity at s = Sc- The name 
double scaling limit |l8l[27lll6lll2lllTl[30l[29l[23l[n] means that we consider a 
regime where the limits s — Sc and N ^ oo are related by a scaling relation 



(s-s,) iV-" = 0(l) 



(1-15) 



where a is some appropriate exponent (a = if Z is not singular). 
We are interested in computing expectation values of resolvents: 



Tr 



xi - M 
as well as in their cumulants 

^ni^Xi, • • • ) Xf-, 



1 ...T. 1 



x„ - M 



X\ Ajj Xfi Aj^ 



Tr 



1 ...Tr ' 



xx - M 



Xn-M 



(1-16) 



:i-17) 



The density correlation functions Pn{xi, ■ ■ ■ , Xn) can be easily deduced from them: den- 
sities are discontinuities of resolvents, and resolvents are Stieljes transforms of densities, 
for example for the 1-point function: 



a)i(a;) 



pi{x')dx' 



X — X 



Pi{x) = 77- {ooi{x - iO) -u;i{x + iO)). (1-18) 

ZtTC 



Imagine, that, for s < Sc, the potential V{x) is such that there is a large expansion 
of the type: 

oo 

lnZ = J2iN/sr'^f„ (1-19) 

9=0 

and similarly: 

oo 

Un{Xi, ...,Xn) = ^(iV/s)2-29-n^(9)(a;^, . . . , ^^). (1.20) 

g=0 

First, let us emphasize that such an expansion does not exist for any potential V, it 
exists only if the integration contour for the Aj's is a "steepest descent contour" for 
the potential V (i.e. a landpath and bridges path in the Riemann-Hilbert language 
of [HI [21] )• For instance it was proved 0^ that such a large N expansion holds for s 
sufficiently small. 
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From now on, let us assume that we are in a situation where such an expansion 
exists. In that case, the coefficients ivn^ and fg were computed in [351 [191 136] , and they 
are the "spectral invariants" of some spectral curve associated to V/s. The spectral 
curve y{x), in that case, is the function y{x) = V'{x)/2 — ujf^\x), it is the large 
density, also called equilibrium density y{x) = iirpeqix) = iTxpf\x): 

Theorem 1.1 (proved in 17^ ) The coefficents fg and Cjn'^ of the topological ex- 
pansion of In Z and Un, o'^e the spectral invariants (in the sense of f3Ef ) of the spectral 
curve: ^ 

y[x) = in pcq{x) = -V'{x) - Ljf\x) . (1-21) 

Remark 1.1 We refer the reader to [36] to see how to compute the spectral invariants of 
an arbitrary plane curve y{x). We shall give an explicit example of computation of spectral 
invariants for formal matrix models below in section fl .2. 11 see theorem 11.31 

Let us say that we shall not be really using any deep result of [36J in this article, except 
the theorem 8.1. of [36j (which is very easy to prove by recursion). 

Here, as is well known in random matrix theory, ujf'\x) is an algebraic curve (hy- 
perelliptical for the 1-matrix model), with typical square-root branchpoints at the 
endpoints of the distribution of eigenvalues, we shall write it: 

cof\x) = ^V'{x)-y{x) (1-22) 

where y{x) is the square root of some polynomial 

y'^ = Polynomial (x). (1-23) 

Generically, this s-dependent polynomial has only simple zeroes, and y{x) has 
square root singularities, but for some appropriate choices of s = Sc, the polynomial 
may have multiple zeroes, and we shall consider that, at s = Sc, there is a zero of order 
2m + 1 at X = 0: 

s = Sc — > y(x)~a;™+i (1-24) 
When s is close to Sc, we typically have: 

m 

y ~ ^x'^^^ (s - s,)^(l + 0{{s - s,)^) ) , (1-25) 

k=0 

which we write: 

m 

y{{s - Sc)~ X) ^ {S - Sc)^ y{x) {1 0{{S - Sc)~) , |/(x) = Cfcx'^+a. 

k=0 

(1-26) 

At s = Sc we have y ~ x"^~^2 and at s 7^ Sc we have y ~ ^/x. 

Notice that a regular endpoint corresponds to m = 0, and in that case, Sc can be 
chosen as any value of s. 
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1.1.1 Double scaling limit 

In this article, we shall be interested in the behavior of tui^^ when s is close to Sc, and 



the Xj's are in the vicinity of a branchpoint. Theorem 8.1. of [52], implies that after 
rescaling, we have (when 2 — 2g — n < 0): 

Theorem 1.2 (theorem 8.1. of J3^): If m > 0: 

f,-{s- (1.27) 

and if m > 0: 

~ (s - 4^)(xi, . . . , x„) (1 + 0((s - s,)^)), (1-28) 

where the fg 's and ujn"^ 's are the spectral invariants of 13^ for the spectral curve y{x) 
appearing in eg. U-26\) . 

Remark 1.2 This theorem is very easy to prove by recursion on n and g from the definitions 
of spectral invariants in [36]. A more detailed proof is also given in section 4.8.2 of [38]. 

Our goal in this article, is to show that the coefficients oon^ and fg can also be 
computed from determinantal formulae for the {p, 2) kernel appearing in the {p, 2) 
minimal model, this is our theorem 12.31 

1.1.2 Specific heat 

If m > 0, we consider the resummation to leading order: 
InZ = J2(N/sr'^fg ~ J](iV/.,)(2-2^)(s - s^f-'^^m^^ fg = F{{s - .JiviSi). 

9 9 

(1-29) 

This shows that the double scaling limit is — oo, s — > Sc with a scaling: 



t= {s-Sc)N^ =0(1). 



:i-30) 



This is a special case of the double scaling limit (s — Sc) ^ N ^P"^"? i)/{p+'?) for general 
We defined the function 



oo 

> 2m + 3 



F(t) = 5^t(-^^)^/,. (1-31) 

9=0 
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Consider its second derivative, often called specific heat: 



X — -\ l-o(2m+3) 

""^^^ = = (1-32) 

9 

(1-33) 

We shall prove in this article, that, as claimed in many physics articles, this function 
satisfies the m+V'^ Gelfand-Dikii non-linear equation. For instance if m = 1, it satisfies 
the Painleve I equation: 

3^2 -]^u = t. (1-34) 

Moreover, it is well known from general considerations in statistical physics, that 
the free energy — In Z should be convex, i.e. u{t) should be negative for t > 0: 

u{t) < 0. (1-35) 



Remark 1.3 The case m = 0, needs some care. The correlation functions a;„'s, indeed 
correspond to the s — > Sc hmits of u)n's, i.e. a zoom x ^ {s — Sc)x near a regular branch 
point, but the free energy F = InZ is not divergent at s = Sc, and thus F{t) cannot be 
seen as the s ^ Sc limit of InZ. In that case m = 0, the 1st Gelfand-Dikii equation is not 
differential, it is simply 

u{t) = (1-36) 

and, if it made sense, it would correspond to a free energy diverging as -P ~ N'^, but in fact 
the free energy F is not divergent, and one finds that fg = for g > 1. 

1.2 Formal matrix models and combinatorics point of view 

As it was discovered by Brezin-Itzykson-Parisi-Zuber [2], matrix integrals are (in the 
formal sense) generating functions for counting discrete surfaces of a given topology 



Consider the potential: 

2 1 f^+l 

V(x) = - 5V(x) , 5V(x) = Y-^' (1-37) 

2s s ^-^ 7 

i=3 ■> 

Formal matrix integrals are defined as: 

Z= [ dMe-TT^^M) 4£,f ^ /■ dM {trSViM))'^ e-^^'"^. (1-38) 

^formal ^ J Hm 
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where dM is the Lebesgue measure on Hn, normahzed such that J dM e'^'^^ ^ = 1. 
In other words, we have exchanged the integral and the summation over k. Z is a formal 
power series in s, where each coefficient is a finite sum of polynomial expectation values 
of a Gaussian integral: 

oo 

Z = l + ^sMj. (1-39) 
After taking the Log, we also have a formal power series: 

oo 

lnZ = ^sMj. (1-40) 

i=i 

It was noticed by t'Hooft [M] and then BIPZ[T^, that, after dividing by A^^, each 
coefficient Aj is a polynomial in l/N"^, namely: 

9max(i) 

g=Q 

Then, one defines: 

oo 

f, = Y.''''"~"' "^3,9 (1-42) 

which is also a formal series in powers of s (one can easily prove that Aj^g = if 
J ' + 2 — 2(7 < 0). In that case, the following large N topological expansion holds as an 
equality between formal series of s, order by order in s: 




(1-43) 

We emphasize that this equality is not a large asymptotic expansion, it is a small s 
asymptotic expansion, and order by order in s the sum over g is finite. 



It was proved by BIPZ in 1978 [14], by a mere application of Wick's theorem, that 
fg is the generating function for maps of genus g: 

„#vcrtices 

f _ \^ „#triangles #quadrangles #(Q!+l)-gons . .n 

maps, genus g 

Similarly one may compute formal expectation values: 

= V -T-T^ (trM'i ... trM'") , , 

Z-^ ^'1+1 ^ ' c,formal 

/ -^1 ■ ■ ■ -^n 
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tr — . . . tr — ) (1-45) 



Xi-M Xn- M 



c, formal 



where c means cumulant, and formal means that we compute the integral by exchanging 
the order of the Taylor expansion of q-^^^ ^^(^^'^ and the integral, as in eg. (11-381) . 
ujn{xi, ■ ■ ■ ,Xn) is thus a formal power series in s, whose coefficients are polynomial 
expectation values of a Gaussian integral. For each power of s, the coefficient in ujn is 
a polynomial in the 1/xj's, and is a polynomial in 1/A^. We write: 

oo 

LO^x,, ...,Xn)= ^(iV/s)2-2^-"c2;i^)(xi, ...,Xn) (1-46) 

9=0 

which is an equality between formal series of s (order by order the sum over g is finite) . 

It was also proved by BIPZ []A\, that oon^ is the generating function for maps of 
genus g, with n marked faces, and with 1 marked edge on each marked face: 

„#vertices „#triangles #quadrangles #{a!+l)~gons 

uj(s) = V - - - • • • (1-47) 

^ #Aut x[^+\..x^^-+^ ^ ' 

maps, genus g in 

where k is the length of the i^^ marked face. 

All this is now standard result in combinatorics of maps. 

1.2.1 Spectral invariants and spectral curve 

More recently, it was proved in ES] , that the functions oon^ and fg are the spectral 
invariants (in the sense of [36]) of the formal matrix model spectral curve informal mm = 
(x, y) defined parametrically by: 

^formalMM - | ^^^^ ^ ^.^^ _ (1-48) 

where the coefficients a,'j,Uk are entirely determined by the following algebraic con- 
straints: 

d d 

x{z) - J2 ^i^y = Yl ^^^^^ + ' ""1 = - ' -"0 = 0, (1-49) 

j=2 k=l ^ 

which give an algebraic equation for a and 7, whose solution we choose such that a 
and 7^ are formal power series of s starting with: 

^^ = s + 0{s^) , a = 0{s). (1-50) 

(We give the example of quadrangulations below) 

Then, from [351 [H] we have: 
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Theorem 1.3 (proved in [7g[ l3B^): a)^f^ and fg are the spectral invariants (in the 
sense of 136^ ) of the formal matrix model spectral curve i^formaiMM = {x,y). 
The spectral invariants of that curve are defined as follows (see fSBjl): 

U^^\x{Zi), X{Z2)) = X-// A - r-^f \ ^ V ^^2 ' (1-^1) 

[Zi - Z2Y x'[Zi)x'[Z2) [x[Zi) -X[Z2)Y 

and with J = {x{zi), . . . we have recursively: 

U}i^ll{J,x{Zn+l)) 

Res ,.;^^^f t , U72\J, i{zlx{z)) 



2x'{Zn+l) ^^±1 Zn+l {x{Zn+l) - x{z))y{z) 
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h=0 ICJ 

(1 - 52) 



and for g > 2: 



1 " -Xg), 



where $'(2;) = y{z) x'{z), and for g = 1: 

1 



/i = ^ln(7V(im-l)) (1-54) 



and for g = 0: 



fo = \{^J2-Yiuj+i-Uj-iy + '^{-iy{u2j-i-u2j+i) 



+ — + sMn ' 



(1-55) 



Example, for quadrangulations we have S4 7^ and all the other Sk = 0. That 
gives: 



{ x{z)=-i{z+\) 



^quadrangulations ^ 



7' = gi^(l- VI-I2SS4), 

uq = U2 = a = 0, 



U3 = -SiY- 



That gives: 

? 1/7^, \2 7'/ \2 / \ 2s7 3s^ 2i 7' 

Jo = 2 (^Y("3 - ^^i) +^(^1) - 2s-f{ui - U3) + —{us) + — + s In 



(1-56) 



(1-57) 
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A = lln(2(2.-7')), (1-58) 

178s=^ - 465sV + 420s7^ - I3O76 
- 6! 4 (72 - 2sr ' ^^'^^^ 

and so on... Notice that at s = Sc = y has a singular branch point y (x — 27)'^/^, 
and when s s^, fg diverges as (s — Sc)^*-^^^^-*. Anticipating on what follows, we see 
that the exponent 5/4 = {p + q)/ {p + q — I) indeed corresponds to the (p, q) = (3, 2) 
minimal model, of central charge c = 0, called pure gravity. In other words, the 
statistics of large quadrangulations is equivalent to the pure gravity (3, 2) conformal 
minimal model, i.e. Liouville field theory. 

1.2.2 Limits of large maps 

It is well known that the asymptotic large size behavior of a number of objects is related 
to the singularities of its generating series. Therefore, the number of maps with a large 
number of vertices (large maps), is governed by the singularities, i.e. the values of Sc, 
such that In Z is not analytical at s = Sc. One should consider the singularity Sc closest 
to the origin, i.e. \sc\ minimal, and see how the fg and Un^'s diverge at s ^ Sc- Thus, 
the scaling limit s — > Sc of a formal matrix integral near a singularity Sc, corresponds 
to the asymptotics of large discrete surfaces. 

For instance, one easily sees from eq. fll-44p . that the expectation value of the number 
of vertices for maps of genus g is: 

< ^vertices >= s In /„ (1-60) 
as 

and thus large maps become dominant when s — Sc a singularity of fg. Typically, if 
we have an algebraic singularity of the type fg ~ (sc — s)""" fg, the expectation value 
of the number of vertices is: 

< #vertices > (1-61) 

Sc - s 

and we see that {sc — s), i.e. the distance to critical point, can be thought of as the 
"mesh size", so that the area (i.e. number of vertices times mesh size) remains finite 
in the limit. 

Another way to say that, is imagine that /„ has an algebraic singularity of type 



f, ~ (Se - s)-°« /, (1-62) 



and notice that 



^ ( 7^ ) i-s/s.r = ± "^^^ is/s.r (1-63) 



(1 - s/sc)-"^ - ^ ^ 
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This means that the (possibly weighted) number of maps of genus g with v vertices, 
behaves for large v as: 

^^^^ v\T{ag) ^''^ T{a,) 

where we used the large v Stirling asymptotic formula for the F-function. 

Similarly, the s — > Sc asymptotics of Cjn\xi, . . . ,a;„) give the enumeration of large 
maps with n marked faces, and if we also rescale Xi — > (sc — s)""'^ Xi, we ca also consider 
large maps with large marked faces. 

Therefore, we see that the enumeration of large maps, is asymptotically given by 
the knowledge of: 

- the exponents an^g (and ag = ao,g); 

- the critical point Sc 

- and the prefactor fg. 

• It turns out that the critical point Sc is independent of n and g, and it can be 
easily found from the resolvent ujf'\ it is not universal, it is strongly model dependent. 

• The exponent turns out to be proportional to {2 — 2g — n): 



a. 



{2-2g~ n) (1 - 7string/2) (1-65) 



where 7string is a universal exponent, it depends only on (p, g), and it is one of the 
exponents computed by the famous KPZ formula. Here, we shall see that it is: 

-2 

7strmg = — -■ (1-66) 

p + q-l 

• The last thing to compute, is the prefactor fg, or ujn'' . Here, in this article, 
we prove in theorem 12.31 the long claim statement that this prefactor is the same as 
the one computed directly from conformal field theory technics, with the Liouville 
theory coupled to matter reprensented by a minimal model (p, q) of central charge 
c = 1 — 6(p — qY /pq (notice that the (3,2) model has c = and thus is called pure 
Liouville gravity). In particular, we show that the generating function of the coefficients 
/g, satisfies the Gelfand-Dikii non linear equation, see eq.( ll-3ip . 

1.2.3 Singularities of spectral invariants 

One can easily convince oneself that the algebraic equations eq.f ll-49p obeyed by a and 
7, are singular whenever y'{l) = or y'{—l) = 0, and then from theorem 11.31 one can 
see that the fg's and cul^^'s diverge whenever y'{l) = or y'{—l) = 0, i.e. whenever y 
doesn't behave as a square root branchpoint. 
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Let us assume that we fix the parameters Sk and s = Sc such that: 

y{z) ~ {x{z) - x{l))"'+^ . (1-67) 
This can be obtained for instance if we choose: 

V'{x) = ix-a-2r (Tm+iix -a)- T„(x - a)) 

., = (-a-2)-(T;+,(-a)-T;(-a)) 

Tm+ii-a) = Trai-a) (1-68) 

where Tm{z + z^^) = z^ + z''"^ is the Tchebychev's polynomial of degree m. In that 
case we have at s = Sc: 

x{z) = z+\ (^.m\ 
y{z) = {z- l)2™+i - (i - l)2-+i 

When s is close to Sc but not exactly equal to Sc, we have like in eq. (ll-26p : 

m 

y((s-Sc)^x) ~ (s -Sc)5^ l/(x) (l + 0((s-Sc)^) , |/(x) = Cfc x'^+i 

fc=0 

(1-70) 

At s = Sc we have y ~ x^^s and at s 7^ Sc we have y ~ a/x — a;(l). The value of m 
and the coefficients Ck depend on which limit of large maps we are interested in. Indeed 
we may fine-tune the coefficients Sj, in order to favor one value of m or another. 
Again, theorem 8.1. of [36] implies that: 

Theorem 1.4 (theorem 8.1. of 13^ ): If m > 0: 

h - {s - s^f-'^^m^^ f,, (1-71) 

and if m > 0: 

'^n\{s - Sc)'^Xi, . . . , (S - Sc)^X„) 

~ (s - sc)('''^-")i^-^ ui^^x,, . . . , x„) (1 + 0((s - Sc)^)), (1-72) 

where the fg 's and ujn^ 's are the spectral invariants of 13^ for the spectral curve y{x) 
appearing in eq.H^TW. 



It was argued and highly debated, that this limit should be equivalent to the Liou- 
ville gravity conformal field theory, coupled to some matter field given by a conformal 
minimal model {p,q) of central charge c = 1 — 6 ^^"^-' . Intuitively, discrete surfaces 
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made of a very large number of small polygons, should give a good approximation of 
smooth Riemann surfaces... 

It was indeed proved that the critical exponents —ag = (2 — 2g) 2^+2 same 
(given by KPZ formula [5T1 [3T1 [32] ) as those of the Liouville conformal field theory, 
but it is only recently that it became possible to compute explicitly partition functions 
and correlation functions on both sides: on the matrix model side (in particular in the 
double scaling limit), and in the conformal theory side. 

On the Liouville conformal theory side, recent progress was obtained following 
Zamolodchikov, Belavin, Hosomichi, Ribault, Teschner, ... [3llH 151 163 1 IGTI H3] . 

On the matrix model side, recent progress was obtained in [35], and formalized as a 
special case of the symplectic invariants of [36] , which allow to compute all correlation 
functions of all genus. 

From here, we can repeat all what was said in section 11.11 after theorem II. 1[ 

In this article, we shall show how to apply the spectral invariants method of [36], for 
the double scaling limit of matrix models which is expected to coincide with Liouville 
theory. 

We prove that the scaling limits of the matrix model correlation functions, i.e. the 
generating functions counting discrete surfaces, is indeed the {p, 2) reduction of KdV 
satisfying string equation, i.e. the minimal model (p, 2). 



2 Minimal models 



There exists several equivalent definitions of minimal models coupled to gravity. They 
correspond to representations of the conformal group in 2 dimensions. They are clas- 
sified by two integers {p, q), and their central charge is: 

c = 1 - 6^^ (2-1) 
pq 

Some of them have received special names: 

• (1, 2) = Airy, c = —2 (related to Tracy- Widom law [66] ) 

• (3, 2) = pure gravity, c = 

• (5, 2) = Lee- Yang edge singularity, c = — ^ 

• (4, 3) = Ising, c = i 

• (6, 5) = Potts-3, c = I 

Minimal models can also be viewed as finite reductions of the Kadamtsev- 
Petviashvili (KP) integrable hierarchy of partial differential equations. 
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The case g = 2 is a little bit simpler to address, and is a reduction of the Korteweg 
de Vries (KdV) hierarchy. 

The KdV hierarchy, and the minimal models {p, 2) have generated a huge amount 
of works, and have been presented in many different (but equivalent) formulations. For 
instance in terms of a string equation for differential operators, in terms of a Lax pair, in 
terms of commuting hamiltonians, in terms of Schrodinger equation, in terms of Hirota 
equations, in terms of isomonodromic systems, in terms of Riemann Hilbert problems, 
in terms of tau functions, in terms of Grasman manifolds, in terms of Yang-Baxter 
equations, ...etc, see [2] for a comprehensive lecture. 

All those formulations are equivalent, and let us recall some of the well known 
features of the (p, 2) reduction of KdV (see [26l [2]), presented in a way convenient for 
our purposes. 

2.1 String equation 

The KdV minimal model {p, 2) with p = 2m + 1 can be formulated in terms of two 
differential operators P, Q of respective orders p and 2, satisfying the string equation: 

[P, = ^ W (2-2) 

Q = (f- 2u(t) , P = dP -pudP-^ + ... , = — 4 (2-3) 

N dt 

is a scaling parameter, which we can send to zero to get the "classical limit". 
The general solution of the string equation eq. (12-21) is of the form: 

m 

P = ^t^.(g^+V2)^ , t^ = l (2-4) 
i=o 

where {Q^^^^'^)+ is the unique differential operator of order 2j -|- 1, such that: 

order[((g^+i/2)^)2 _ g2,+i] < 2j. (2-5) 

For example: 

{Q"')+ = d , (q3/2)^^^3_3^^_^^ ^2-6) 

(g^/')+ = rf^ - bud"" -—d^-—d-^d-—u —. (2-7) 

It is a classical result (see [27]) that it satisfies: 

[{Q'-"')~.M = ]^j^{RMm (2-8) 
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where the right hand side is a function (a differential operator of order 0), and the 
coefficients Rj{u) are the Gelfand-Dikii differential polynomials [27]. They can be 
obtained by the recursion: 



Rq — 2 

The first few of them are: 



R:^ 



Rj+i = —2uRj — iiRj + 



1 



Ri- 



Ri = -2u 
R2 = ^u^ - 



2 7V2 



U 



2iV2 



uu + 



4iV2 



u 



u 



(2-9) 



(2-10) 



and in general: 



R,{u) 



2 {-ly (2j - 1)!! 



u 



(2i~2) 



(2-11) 



J! (2iV)2i-2 
After substitution of eg. (12-41) into the string equation eg. (12-21) . the property eg. (12-81) 
gives a non-linear differential eguation for the function u{t): 



^tjRj+i{u) =t. 

j=0 



(2-12) 



Since Rq = 2, we see that we can identify t with t = — 2t_i. 

• For instance for Airy p = 1, this gives: 

-2u = t. 

• For instance for pure gravity p = 3, this is the Painleve I eguation: 

1 



3 m" 



2N^ 



u — 2tQU = t. 



For instance for Lee- Yang p = 5, we have: 



-uii -ii -u+ti{3u -u)—2tr)U = t. 

2m 4A^2 8]V4 ^ IV 2A^2 ^ " 



(2-13) 
(2-14) 
(2-15) 



2.2 Tau function 



We define the Tau-function r(t, to, • • • , ^m) and its log, the free energy function 
F{t, to,. . ., tm) = In r(t, to,. . . , tm) such that: 

N-^F = u. (2-16) 

The Tau-function has many other properties, which can be found in textbooks and 
classical works on the subject [21 [551 [56], but which are beyond the scope of the present 
article. 
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2.3 Lax pair 

Consider the following matrices: 



and for any integer k: 



Vk{x,t) 



1 

X + 2u{t) 



Ck —Ak 



(2-17) 



(2-18) 



where Ak{x,t), Bk{x,t),Ck{x^t) are polynomials of respective deg 
which are determined by: 



Bk{x,t) = Y,^''~' R3 

j=0 



Ak 



2N 



Ck = {x + 2u) Bk + ^Ak. 



The recursion relation eg. (12-91) implies that Bk satisfies the equation: 
2uBk + 2{x + 2u)Bk 



^ Bk = -2Rk+i{u) 



2A^2 



and we see that the matrix Vk{x,t) satisfies: 



l_d_ 

Ndi 



the right hand side is independent of x, and is proportional to ■^TZ{x,t). 
2.4 Lax equation 

If we consider u solution of the string equation eq. (l2-12p . then, the matrix: 

m 

V{x,t) = ^tj'Dj{x,t) , tm = 1 

j=0 



satisfies the Lax equation: 
1 d 



2 d 



^ Qt^M + [^^(x,^),7^(a;,^)] = --— 7^(x,t) 

which can also be written as a commutation relation: 



This relation means that the operator ^ + 'D{x,t) is a Lax operator [2J. 



(2-19) 



(2-20) 



Vk{x,t) + [Vk{x,t),n{x,t)] = -^Rk+i{u) (^^ (2-21) 



(2-22) 



(2-23) 



(2-24) 
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2.5 The differential system 

The Lax equation eq. fl2-24p is the compatibihty condition, which says that the following 
two differential systems have a common solution \E'(x,t): 

l^-^^^J(x,t) = -^Vix,t)^ix,t) , 1 A^(a;,t) =7^(a;,t)^(x,t) (2-25) 
and 'if{x,t) is a matrix such that: 

= (t , det^ = 1. (2-26) 



In particular we have the Schrodinger equation for ip: 

^{x, t) = ix + 2uit)) ^(x, t) (2-27) 

where t can be interpreted as the space variable, and x the energy, x is called the 
spectral parameter. 

2.6 Correlators 

Consider the Christoffel-Darboux kernel associated to the system V{x): 

K{x„x,) = (2-28) 

Xi — X2 

Definition 2.1 We define the connected correlation functions by the " determinantal 
formulae" : 

Wi{x) = lim K{x, x') ^ = ^'(x)0(x) - 'ip'{x)(p{x) (2-29) 

and for n >2: 

Wr,{xu...,Xn) = -—^^-{-ir n^(^-W'^-(^+i)) (2-30) 

\Xi X2) . 

For example: 

W3(xi, X2, X3) = K{xi, X2)K{X2, X3)K{X3, Xi) + K{xi, X3)K{x3, X2)K{x2, Xi). (2-31) 

Although we have not written it explicitly, the kernel K and the correlators Wn 
depend on t. 

The non-connected correlation functions are defined by: 

k 

Wn,n.c.{^l,...,Xn) = J2 Yl n^l-^^K-^^)' (2-32) 

k JiUJ2U...UJfc=J «=1 
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where J = {xi, . . . ,Xn} and the sum runs over all partitions of J into k non-empty 
disjoint subsets. In other words, the connected WnS are the cumulants of the non- 
connected ones. 
For instance: 

W2,nAxuX2) = W2iXi,X2) + W,iXi)Wiix2), (2-33) 



M^3,n.c.(a;i,X2,X3) = W3{Xi,X2,X3) + Wi{Xi)W2{x2, X3) + Wi{x2)W2{Xi, X3) 

+Wiixs)W2ix,,X2)+Wiixi)Wiix2)Wi{xs). (2-34) 

The formula eq. fl2-3Up is called "determinantal formula" , because for the non-connected 
correlation functions we have: 

Wn,n.c.{^ij • • • , a;„) = det(i^(xi, xj)), (2-35) 

where det' means that when we compute the determinant as a sum over permutations 
of products (—1)'^ K(xj, Xo-(i)), then if a{i) = i we replace K{xi,Xi) by Wi{xi), and 
if a{i) = j and a{j) = i, we replace K{xi,Xj)K{xj,Xi) by —W2{xi,Xj), see [7]. 
For instance W^^n.c. is the sum of 6 terms coming from the 6 permutations: 

(K{xi,Xi) K{xi,X2) K{xi,X3)\ 
K{x2,Xi) K{X2,X2) K{X2,X3) 
K{x3,Xi) K{X3,X2) K{X3,X3) J 
= Wi{xi)Wi{x2)Wi{x3) +Wi{xi)W2{x2,X3) + Wi{x2)W2{xi,X3) 
+ Wi{x3)W2{xi,X2) + K{XI,X2)K{X2, X3)K{x3, Xi) 

+K{xi,X3)K{x3,X2)K{x2,Xi) (2-36) 



It was proved in [8j, that the correlators Wn satisfy an infinite set of equations, 
called loop equations, and equivalent to Virasoro constraints for the r function. The 
loop equation simply states that the following quantity: 

Theorem 2.1 Loop equations (proved in JE/)-' 

PniS^i -^li ■ ■ ■ 1 ■^n) 

n.c. (•^5 -^1 ) • • • 5 •^n) 
_|_ d Wn{x, Xi, . . . , Xj-i, . . . , Xn) — Wn{Xi, . . . , X„) 

Z—/ f)nr- ■ 



dXi X — X 

(2 - 37) 



j = l 1 



is a polynomial of the variable x. 
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For example, one can easily check that: 

Pq(x) = W2{x, x) + Wi{xf = - det V{x, t). (2-38) 

Notice that 

Id 

--detV{x,t)=2Bix,t) = 2Y,hB,{x,t). (2-39) 

i=o 

2.7 Example: Airy kernel 

Let us write the (1, 2) model, i.e. m = 0. We have: 

P = d , Q = d?-2u (2-40) 

the string equation is: 

1^.01 = -^« = ^ (2-41) 

i.e. 

u{t) = -^-= t_i (2-42) 

The Lax pair is: 

The differential system is: 

1 

1^ dx^^^'"' \x-t 



^ vl>(x,t) = -f ^ 1 ) vl;(a:,t) (2-44) 



I.e. 

r2 



^" = N^[x - t)%lj (2-45) 
whose solution is the Airy function [1]: 

i){x,t) = Ai{Ni{x-t)) , Vi(x,t) = -Ai'(ivi(x-t)) (2-46) 
and the other independent solution is the "Bairy" function [Ij: 

(l){x,t) = Bi{Ni{x-t)) , ^{x,t) = -Bi\Ni{x-t)) (2-47) 
and thus the kernel is the famous Airy kernel |65j : 



K^.y{t + N~^/^x„t + N~y^x,) = -4^XxOE.(x.) - A^ix^)B^'ix,) ^^-48) 

Xi — X2 

The Airy kernel plays a very important role in many problems, in particular in the 
universal laws of extreme values, related to the Tracy- Widom law [66] . 
The r function is simply: 

T = e- — . (2-49) 
For the Airy system, the polynomial of theorem 12.11 is simply: 

Pn{x) = 4(x + 2u) 6n,o. (2-50) 
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2.8 Classical limit 



The classical limit is the large N limit, or equivalently, it is also the large t limit. 

Intuitively, in the classical limit, P and Q commute, and they can be represented 
without differential operators. In this limit d —>■ z can be represented as a number, 
and operators Q = d? — 2u and P are replaced by functions of z and t. Therefore, in 
analogy with Q = — 2u{t), and P = d"^ + . . ., let us define two functions x{z, t) and 
y{z,t): 

x{z, t) = z^ - 2uo{t) , y{z, t) = zP + ... . (2-51) 

In the classical limit, we replace the string equation [P, Q] = N^^ with a Poisson 
bracket: 

^y^""^-^- dz dt dtdz ^^^^^ 

whose general solution is: 

x{z, t) = z'~ 2u,{t) , y{z, t) = Y^ t, h^^' (1 - ^).•+V2^ ^ (2.53) 

i=o ^ ^ ^ + 

where ()+ means the positive part in the large z Laurent series expansion. Explicitly 
we get: 

i=0 1=0 J \ J ' J 

The string equation {y, x} = 1 reduces to: 

wol/'(0) = ^, (2-55) 

i.e. 

f^t>o(-W2)^-^^^ = -^ (2-56) 
which can be integrated with respect to t and gives a polynomial equation for uo{t): 

V(u,) = ±t,(-u,/2y*^ ff±^^^l (2-57) 

which is clearly the classical limit of eq. (12-121) . In other words, the non-linear differential 
equation eq. (l2-12p for ■u(t), becomes an algebraic equation for uo(t). 

For example, for pure gravity m = 1 we have the classical limit of eq. (12-141) : 

AV{uo) = 3ul- 2to uo = t. (2-58) 
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2.9 Topological expansion 

We now have the polynomial equation eq. fl2-57p : 

V{uo) = t/A (2-59) 

which implies: 

and in general, any derivative of Uq with respect to t can be written as a rational 
function of Uq- 

Since 'Uo('^) satisfies the string equation eq.f l2-12p at = oo, the full solution u(t) 
to the string equation eq. fl2-12p . can be expanded as an A^"^ power series: 



u 



(t)=Mo + 5^iV-2fcMfc(t) (2-61) 



where all coefficients Uk are rational functions of Mq (their denominator is a power of 
V'iuo)). 

For example for pure gravity m = 1, the Painleve equation eq. (12-141) implies that 
to the first few orders we have: 

«W = "0 - ]^ ^6^0 - 2to)-' + 0{N-^). (2-62) 

And the Free energy F{t) such that u = ^ F, also has a 1/A^^ expansion: 

oo 

\nT = F = Y,N^~^'Fg{uo) , Fg=ug. (2-63) 

5=0 

Also, since the coefficients of the differential system V{x,t) depend on u(t), the 
matrix V{x,t) has a 1/A^^ expansion: 

V{x, t) = J2 N-^W^^^ (x, t) (2-64) 

9 

To leading order we have: 

m-,t)=(. ^ ^, ^ (2-65) 
\{x + 2uo) B{x,uo) / 



Bix, «o) = 2 f; t,x^-' 4 ^ '^'^'f (2-66) 



j=0 fc=0 
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Notice that: 

zB{z^ ~2uo,uo) =y{z,t). (2-67) 

The classical spectral curve is given by the eigenvalues of T)^^\x,t)), i.e. the values 
of y such that det {y — T)^'^\x,t)) = 0, i.e., if we parametrize x as x = z"^ — 2uo, we 
have: 

y = ±yiz,t) (2-68) 

where y{z, t) is the function defined in eg. fl2-54p . This explains why we call the function 

y{z,t) the classical spectral curve. 

Written in a parametric form where Uq = Uo{t), the classical spectral curve is thus: 



x{z) = z^ — 2uo 



(2-69) 



It is important to notice that it is a genus hyperelliptical curve, which is equivalent 
to saying that it can be parametrized by a complex variable z (higher genus would be 
parametrized by a variable z living on a Riemann surface), and which is equivalent to 
saying that the polynomial y^, written as a polynomial in x, has only one simple zero, 
located at x = —2uq, all the other zeroes are double zeroes. 

2.10 BKW expansion 

Similarly, we can look for a BKW asymptotic solution of the solutions il){x,t) of the 
differential system. It takes the form: 

^(x, t) ~ — 1 + E «o) (2-70) 

V2(-x-2uo)4 \ k J 

Vi(x,t) ~ e'^^-2-o^''" (x + 2uo)2 ||l + ^Ar-'^Vifc(x,Mo)j (2-71) 

and we recall that z = {x + 2uq) 5 . the BKW expansion of the other solutions and 
0, are obtained by changing the sign of the square root z —z. 
We have the following Lemma: 

Lemma 2.1 Each ipk{x,uo) and ■ipk{x,uo) is a polynomial ofl/z. 

proof: 
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The proof uses the Schrodinger equation eg. (12-271) : 

1 



^lj{x,t) = {x + 2u{t))ip{x,t). 



Let us write: 



/ r* 

ip{x,t) = V/(a;,t)e-' /(-•*'). 

The Schrodinger equation imphes that: 



N\x + 2u{t)) fix, t) = ^fix, t)fix, t) - ^/(x, ty + 1, 



(2-72) 
(2-73) 

(2-74) 



and after differentiating once more with respect to t, we obtain a third order hnear 
equation for /: 



{x + 2u{t))f{x,t)+u{t)f{x,t) 



2 m 



f{x,t). 



(2-75) 



To leading order we have u{t) = uo{t), and recall that u{t) has a expansion, 
therefore, one easily sees that: 



f{x,t) 



N^x + 2uo{t) 



(2-76) 



and by an easy recursion, we see that each fk{x,t) is a polynomial in I/2; with z 
^x + 2uo{t). 

Then, notice that the Poisson equation eq. (12-520 implies: 



dt 



And therefore: 



This implies that: 



x'{z) 2z 2 + 2mo 
d ydx 



dt 



-z. 



f{x,t) 



N / ydx + J2N'-''9k{x,t) 



(2-77) 



(2-78) 



(2-79) 



and where all coefficients gk{x,t) are polynomials of 1/z. 

Since "^{x, t) = v/(x, t) e , we find that ip{x, t) is of the form: 



'ilj{x,t) 



[x + 2uo)* 

where each ijjk{x,t) is a polynomial in 1/z. 



r + Mx,uo) 

k 



(2-80) 
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The proof for tp{x,t) works in a similar manner. □ 

This lemma implies that the kernel also have a expansion: 

N J^^i ydx 

K{z,,Z2) = ^ ^] (l + ViV-^ir,(^i,Z2)), (2-81) 

where each /^^(zi, Z'2) is a polynomial in Xj Z\ and in 1/-22- 

This implies that the correlators also have a expansion: 

Lemma 2.2 

. . . , x„) = ^ iv2-2^?-" . . . , x„) (2-82) 



where each Wn^ is a rational function of the Zi = y/xi + 2uo, with poles only at Zi = 0, 

Wi"^ =y{z,t) (2-83) 



except W2^^ and W^^^ which are: 



Ty{0) ^ _J: I 1 ^ 1 (2 QA) 

' Az,Z2{z,-Z2y {zl-zlf AzMz, + Z2y ^ ^ 

The important point, is that each wlf^ has no other pole than = 0, in particular, 
has no pole at the other zeroes of y{z,t). 
proof: 

Notice that in the products Y[iK{za(i),a(i+i)), all the exponentials cancel, and the 
result is, order by order in N~^, a rational fraction of the Zj's having poles at Zi = 0, 
or at Zi = Zj. Except for and W2^\ the poles at Zi = zj are simple poles, and it 
is easy to see that in the sum over permutations, the residues cancel, therefore, each 
is a rational function of the ^j's having poles only at Zi = 0. The cases of W2 and 
Wi need to be treated separately, and are easy. 

The fact that Wn has a 1/ A^^ expansion instead of 1/N comes from a simple sym- 
metry argument. In the expression of Wn, changing —>■ (p and —>■ (p, can also be 
obtained as changing the order of the Xj's, and since we take a symmetric sum, only 
the terms which are invariant under the exchange ip —>■ (p and ip ^ (p contribute to 
Wn- Exchanging the two solutions ip ^ (p and ip ^ (p, is also equivalent to changing 
—N, and therefore Wn has a given parity in A^. □ 

2.11 Symplectic invariants 

It was found in 18], that the correlators obtained from the determinantal formulae 
eq. fl2-29p . eq. fl2-30p of a Christoffel-Darboux kernel K of type eq. (12-281) . do satisfy loop 
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equations, i.e. for any n and g, and J = {xi, . . . , Xn}, the following quantity: 



h=0 ICJ 



^ d Wi,'\x,j/{X,})~WI,'\X„J/{X,}) 

^ dXn 



=1 ^-"J 



(2-85) 



is a polynomial in x. This property, as was proved in [5], is a direct consequence of 
eq.dMSl) and eq.dMHD, eq.dESO]). 

Moreover we know from section [2.10^ that Wn\x{zi)^ . . . the following dif- 

ferential form: 

W(^)(zi, . . . , z„) = ■ ■ ■ , ^(^n)) ^'(^i) • • • x'{z^) + (2-86) 

is a symmetric rational function of all its variables, and if 2(7 + n — 2 > 0, due to lemma 



12. 2[ it has poles only at = 0, and 



1 



(2-87) 



It was found in [35], that the unique solution of loop equations eq.f l2-85p which has 



a topological expansion for which the wi^^'s have the poles given by lemma [2^21 can 
be obtained by the following recursion relation: 

Theorem 2.2 



>vi^i(zi,...,z„, z„+i) = — ^ Res 



dz 



4 {zl^^ - z^) y{z) 



Wi%-'\z,-z,J) 



a I 



+ E E ^lU'^ ^)<;-Vi(-^' Ji^) (2-88 



/i=0 IdJ 

where J = {zi, . . . ,Zn}, and J2hYl'i' i^^ans that we exclude the terms {h,I) = (0,0) 
and {h, I) = {g, J). 

proof: 

The proof proceeds exactly like in [35]. Write the Cauchy residue formula: 



5 Zn: Zyi-\-1 ) 



dz 

Res 

z^z„+i Z — Zn+1 



Wi'Uzu---,Zn,z) 



and move the integration contour, to enclose all the other poles, i.e. only z 
thus: 

dz 



(2-89) 
0, and 



)^nll{Zl, . . . ,Zn, Zn+l) 



Res 

z^O Zn+l - Z 
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>vSi(zi,...,z.,z) 



Res ^^M^ Wl,i,{x{z,),...,x{z^),x{z)) 

Zn+1 - Z 

(2 - 90) 

Then, insert in the right hand side eq. fl2-85p : 
-2<)(x)iyii(a:i,...,x.,x) = j^j^w[%^{xj)w[l-\{x, J / 1) 

h=0 ICJ 

d w!,'\x,J/{x,})-W!,'\x,,J/{x,}) 



T — 

^ dx^ 



-P^^\x-x^,...,Xn) (2-91) 

and notice that the polynomial Pn'' has no pole and doesn't contribute to the 
residue. All what remains is eq.( l2-88p . □ 

The recursion relation eg. fl2-88p is precisely the definition of the symplectic invari- 
ant's correlators defined in [36]. In [36|, it is explained how to associate an infinite 
family of >Vn^'''s, to any spectral curve defined by a pair of functions {x{z)^y{z)). 

Examples: 
• eg. fl2-88p gives: 



Wf\zi,Z2,Z,) 



eg. flMHj) gives: 



W^^\zi) = -p- Res 
-1 



1 


1 


2y'{0) z\ 




dz 








( - 





(2-92) 



32 dz^ \{zl-z^) y{z) 
y"'{0) 1 1 1 

48?/'(0)2 ^ ~ 16 y'{0) 4 

2.12 Double scaling limit and {p,2) kernel 



2 = 



(2-93) 



By comparison with theorem II. 2[ we conclude that: 

Theorem 2.3 the s ^ Sc double scaling limit of (possibly formal) matrix integrals 
correlation functions w!f^ are the determinantal formula correlation functions of the 
{p, 2) kernel: 



Wi^\x„ 



(2-94) 



where 



a;<f)(xi, ...,Xn) = lim(s - '^)lZXl + i2;(9) (^(^g _ . . . , (s - sj^+ia;^) 

S^Sc 

(2-95) 
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is the double scaling limit of matrix integrals correlators, and w!f\xi, . . . , x„) is the g^^ 
term in the BKW expansion of the determinantal correlator of the (2m + 1, 2) minimal 
model fdef\2l\). 

And similarly ifm>0: 

lim {s - s.Y^'-'^wfg = fg = Fg (2-96) 

where fg are the free energies of the matrix model and Fg are the free energies of the 
(2m + 1, 2) minimal model: 

luT = ^N^-^^Fg. (2-97) 

9 

The double scaling limit is: 

N~^oo , s^Sc , {s-Sc)N^ =0(1). (2-98) 

Therefore, we have proved that, as announced, the double scahng hmit of matrix 
models, is given by the Liouville minimal models (2m -|- 1, 2) coupled to gravity. 

2.13 Parametrics of orthogonal polynomials and Baker- 
Akhiezer functions 

Many approaches of matrix models use some orthogonal polynomials (see Mehta [57]): 

Pnix) = X" + . . . , j Pnix) Pmix) e'^ ^^^^^ dx = hn Sn,m, (2-99) 

which we prefer to make orthonormal: 

ijnix) = —^=- Pn{x), (2-100) 



as well as their Hilbert transforms: 



0„(x) =e^sW I — dx. 

The matrix 



X — X 



(2-101) 



m^{x)=[ ) (2-102) 



i)n-l{x) 0„_l(x) 

satisfies a 2 x 2 differential system Vn{x) with polynomial coefficients of degree at most 
degV^': 

^-^^n{x) = Vn{x) ^„(x). (2-103) 

A* ax 

It is easy to see that tr Vn{x) = and det "^nix) = ^/hn-i/hn is constant. 
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Since (pnix) is discontinuous across the integration contour of dx', the matrix \l/n(a^) 
is also discontinuous, and has jumps of the form: 

vI/„4x) = M/„_.(a;) (J '^Y)- (2-104) 

Therefore "^nix) satisfies an isomonodromic Riemann-Hilbert problem (the jump ma- 
trix, called the monodromy, is independent of n and of V{x)). 

A general method was invented [2l| to find large N asymptotic solution of isomon- 
odromic Riemann-Hilbert problems. In the case where we approach a (2m + 1, 2) 
singularity, the method of Deift& co [23] requires to have an ansatz for a parametrix 
asymptotics for \E'n(a;) in the vicinity of the singularity. 

We claim that the correct parametrix for the (2m + 1, 2) singularity, is the matrix 
of Baker- Akhiezer functions of eg. (12-261) for the (2m + 1, 2) minimal model: 

^n((s - Sc)^x) ~ ^{x) (1 + 0{{s - Sc)^)). (2-105) 

This should be checked by the steepest descent Riemann-Hilbert method of [21]. 



3 Kontsevich's integral 

In this section, we also propose a combinatorical interpretation of the coefficients of 
the fg's and oon^'s, based on the comparison with Kontsevich integral. 

The spectral curve eg. (12-691) of the (2m -|- 1, 2) minimal model 

{x{z) = - 2Uq 
yi^) = t,z'^^'-'' {-uo/2y mpi (3-1) 

is of the same form as the Kontsevich integral's spectral curve (see in [37|, or see 
below), and thus it has the same correlators and spectral invariants Fg as those of 
the Kontsevich integral. Therefore, the correlators and F^'s of the minimal model 
(2m + 1,2) can be written as integrals of tautological classes on the moduli spaces of 
Riemann surfaces. This can be viewed as a proof that the double scaling limits of 
matrix models, i.e. the limit of large maps generating function, indeed coincides with 
topological gravity, as claimed by Witten [7T1 [28] and then proved by Kontsevich [52] . 

The Kontsevich integral [52]: 

Zk{K) = j f/Me-^T^ = e^. , n = ^Tr A"^ (3-2) 
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(to simplify we assume ri = here) is the generating function for intersection numbers 
of cotangent hne bundles at marked points of Riemann surfaces of genus g: 



(g) 



di 



ipi = Ci(A) (3-3) 



M„ 



where ipi is the Chern class of the cotangent line bundle at point i, and where 

n = (21 - ly.l T2i+i (3-4) 

and more generally, correlation functions of the Kontsevich integral give access [37] to 
integrals of Mumford k characteristic classes |60j : 



do 



do+di+.. 

2di + 1! dz; 



(3-5) 



n 

i=l 



do+di+...+d„=dg^n k=l bi + ...+bf,=do,bi>0 
k k n 

di 



1=1 1=1 i=l 



The class kq is the Euler class, and 27r/ti is the curvature form of the Weil-Petersson 
symplectic metrics. The dual times fk are closely related to the r^'s, see the relation 
in [37] or eg. (I3-9P below. 

It was shown in [301 EZ], that Kontsevich's integral's W^^^ = F^'s and correlators 
can be computed as the symplectic invariants Fg = FglSx) of a spectral curve: 



(3-6) 




We see that the minimal model i^(p,2) spectral curve eq. (12-691) can be identified with 
Kontsevich integral's spectral curve Sk, under the identification of times: 

Sk,o-^r,,^3-j^^:^X^U^k{-uo/2) -jr^i^^ (3-7) 
In particular 

1 - ^ rs = y'iO) = (3-8) 
2 —2 Mo 

The dual times are given by their generating function g{r) = f^r^ and g{r) = 
— In (1 — g{r)) with: 

1 - g{r) = -2uo E E ^'^^ (-0^' ^^\uittv^' = ^"'^^^ ^^-9) 
fe>0 I + 
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^W = $^^fcr'= (3-10) 



fc 



I.e. 



j 1=0 •> 



j •'' 1=0 

(3-11) 

4 Derivatives 

The general method to compute derivatives of Fg and Wn^^^s with respect to any 
parameter entering the spectral curve is explained in |36j . 

Here, our spectral curve i^(p,2) depends on the parameters tj's and t = — 2t_i. [36] 
says that we first have to study the variation of y{z)x'{z) under variation of any such 
parameter, and we write it: 

^.'(.)-^,'W = a;w (4-1) 

Here, we find for j > 0: 



and for j = —1: 

z 



m 



A_i(^) = -2uoy{z) = J2t,A'^{z) (4-3) 

j=0 

The theorem 5.1 of [36], then shows that those functions are such that: 



dWi'\z^, ...,_Zn)_^^^^ Wii.iz,, . . . , z„, z,+i) A,-(z„+0 x'(;.„+i) (4-4) 



and in particular for n = 0: 

" Res W['\z) Kj{z) x'{z) = Res w['\z) Aj{z) (4-5) 



and: 



"^'^ = Res . . . Res Wt\z,, ...,Zk) A,,{z,) . . . A,,{zk) (4-6) 
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4.0.1 Example 

The recursion relation eq. fl2-88p gives: 



1 1 



2y'{0) z\zlzl 



(4-7) 



This implies: 



dtj-^ dtj2^ dijj:^ 



Res Res Res }V^^\zi, Z2, z^) A j^{zi) A j2{z2) A j^{zs) 

oo oo oo 

-1 



(4-8) 



^a;,(o)a;,(o)a;3(o) 



Notice that 



A'(0) 



(2j + l)! 



' j!(j + l)!^ 
In particular this implies that: 



-2j 



A-i(O) = 1 



dt^ 2y'{0) 



Uo{t) 



(4-9) 



(4-10) 



and thus, as expected we recover: 




(4-11) 



4.0.2 Example Wl 



(1) 



The recursion relation eq. fl2-88p gives: 



Wr{z) = W^'>{z)x'{z) 



1 



Kr3-2) \z^ {r,-2)z^ 
-iiQ I 1 iio 

z 



4 3^2 E^'+i(-"o/2) 



(4 - 12) 



I 



(2/ + 3)! 
/!(/ + l)! 



and thus: 



Res y^^^\z)AAz) 

oo 

31^ (a;"(0)-2.oA;.(0) E^'^^(-o/2)'|^ 



(4 - 13) 
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In particular: 



~dt 



(2/ + 3)! 
/!(/ + !)! 



(4 - 14) 



4.0.3 Other examples 



r5! 

2! 



{z\ < > +Zl < ri >) + ^rzfzi < iIj,^P2 > 



+T1Z1Z2 < Kitpi > +fiZ^Z2 < Kli'2 > + 2 ''"1 -2^1 -2^2 < /«1 > 



+T2Z^Z2 < K2> 



Kra - 2Yzlzl 



(rg - 2)2(54 + ^zt + ?,zlzl) + QtIzIzI 



(ts - 2) (6X52:^^2 + Qnztz2 + ^T-jZ^Z^^ 



(4-15) 



Wf\z 



(4 - 16) 



1 



128(2 - T,y^^ + I2rlz\2 - r,){50rrz' + 21 r,) 

+z\2 - r3)2(252 + 348 t^tjz^ + 145 r^z'^ + 308 nTgz^) 
+z^{2 - r3)(203 rg + 145 z^r^ + 105 z\ + 105 ;zVii) 
+105 (2 - ra)^^ 



(^1,^2,^3,^4) = 12 



and so on ... 



4.0.4 Homogeneity relation 

Theorem 4.7 of gives another relation which we can apply here: the homogeneity 
equation. Let $(z) such that 



<^\z) = yiz)x\z) 



A-i 

-2un 



2z = Y,t,K] 



j=0 



and thus: 



3=0 



(4-18) 



(4-19) 
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We have: 



(2 - 2(7 - n) W(^)(zi, . . . , z„) = Res . . . , z„, $ 

(4 - 20) 

In other words, is homogeneous of degree 2 — 2g — n. 

5 The (2m, 1) minimal model 

Another kind of universal hmit of matrix models may arise when two connected compo- 
nents of the eigenvalues support merge, typically the equihbrium density of eigenvalues 
behaves as: 

y ~ x^"". (5-1) 

The case m = 1 was treated in [TOl [111 [13 The results concerning general m, were 
described without proof in [12]. The universal limit is given by the (2m, 1) reduction 
of the mKdV hierarchy. All the results can be proven in a way very similar to the 
(2m + 1, 2) case, and here we just summarize the results the results stated in [T2] . 




Figure 2: When two cuts merge, the equilibrium density of eigenvalues behaves like 
a;^*" near the merging endpoints. The universal eigenvalues statistics is given by deter- 
minants of the (2m, 1) kernel, associated to the integrable mKdV hierarchy. 



36 



5.1 The (2m, 1) minimal model 

Again, we shall write a Lax pair T>{x,t) and TZ{x,t) where T>{x,t) is polynomial in x 
of degree 2m, and TZ{x,t) is polynomial of degree 1, and such that: 



= 



We choose: 



n{x,t) 



x + u{t) 
-x + u(t) 



and the matrix T)(x, t) is of the form: 

V{x,t)^J2^k'Dk{x,t) 



(5-2) 



(5-3) 



(5-4) 



(5-5) 



with: 

. .^ _ f -Ak{x,t) xBk{x,t) + Ck{x,t) 

'^'^^^''^ \xBk{x,t)-Ck{x,t) Au{x,t) 

where A^, 5^, Ck are even polynomials of of degree degA^ = 2k — 2, deg Bj- = 2k — 2, 
deg Ck = 2k. They can be found by recursion: 



Ck+i 
Bk+i 

A-k+i 



, So = 
= x^Ck + Rk{u) 
= x^Bk + Rk{u) 



Co = l, 



(5-6) 

where Rk{u) and Rk{u) are the modified Gelfand-Dikii differential polynomials: 



^o(«) 

-Rfe+i(^i) 
Rkiu) 



u , -Ro(w) = 
uRk{u) - ^Rk{u) 
uRk{u) 



u 



(5-7) 



For example: 



^i(«) = y-| , ^i(^) 
Mu) = ^ ^ , R2{u) 



3u^ 




~8~ 


T 








8 










+ 32 





bu ii 



uu 



(4) 



16 



16 



uu 



(5-8) 
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The matrix T>{x,t) satisfies 

[d^-v{x,t),dt-n{x,t)] = o 

if and only if u(t) satisfies the string equation: 



(5-9) 



'^tkRk{u) = -tu. 



k=0 



The Baker- Akhiezer functions 

^(x,t) 



f{x,t) (f){x,t) 
^/j{x,t) (f){x,t) 



are given by the common solutions of the two compatible systems: 



^^^(x,t) =P(x,t)^(x,t) 



^^^(x,t) =7^(x,^)^(x,^). 



(5-10) 



(5-11) 



(5-12) 



It was claimed in [12j that the parametrics asymptotics of orthogonal polynomials near 
the singularity are given by: 



COS ((n -|- 2)^^) ~ sin ((^ + 2)''^^^ 

1 



j_u{t) cosvre / 

'^ ^(sinTre)^ ( cos (3(n + -)vr6) ^(x) 



(5 - 13) 



— sin (3(n H — )7re) ^ipix) ) -|- . . . 



and where e/(l — e) is the ratio of the number of eigenvales in the 2 cuts which merge 
at the singularity (e = 1/2 is the symmetric case). 



The Christoffel-Darboux kernel K{xi^X2) is the same as eq. fl2-28p : 

V'(Xi)0(x2) - ^(Xi)</)(X2) 
K{Xi,X2) = 



(5-14) 



Xi - X2 

and the correlators are obtained by the same determinantal formulae eq. (12-301) . 

And again the claim is that the determinantal correlators of the (2m, 1) minimal 
model, are the limits of matrix models correlators. 



6 Conclusion 

In this article, we have summarized some properties of scaling limits of matrix models 
(formal or not), known for a long time. We have provided a mathematical proof that the 
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asymptotics of the a)„ s and /^'s are indeed those obtained from conformal field theory. 
Our proof is based on the fact that the hmits ujn^^s of matrix models correlators, are 
the spectral invariants of the limit spectral curve, and the fact that the determinantal 
correlators of the (2m + 1,2) minimal model kernel are also the spectral invariants of 
the same spectral curve. 

We recall that those correlators can be interpreted in the Kontsevich integral's 
framework, and have a combinatorial interpretation as intersection numbers of some 
tautological classes on the moduli spaces of Riemann surfaces. 

We claim that the same method can be applied to other sorts of universal limits, 
in particular the merging of two cuts like in [12], and hopefully, we can work out the 
same kind of proof for multi-matrix models, whose universal limits should be the {p, q) 
minimal models with arbitrary p and q. Unfortunately, one of the key points should be 
the equivalent of [S] (i.e. the fact that determinantal correlators obey loop equations), 
but this is not proved yet for differential systems of order q > 2. 
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